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Abstract 

Water soluble quinones are a group of cytotoxic anti-bacterial compounds that are secreted by many species of plants, 
invertebrates, fungi and bacteria. Studies in a number of species have shown the importance of quinones in response to 
pathogenic bacteria of the genus Pseudomonas. Two electron reduction is an important mechanism of quinone 
detoxification as it generates the less toxic quinol. In most organisms this reaction is carried out by a group of flavoenzymes 
known as NAD(P)H quinone oxidoreductases. Azoreductases have previously been separate from this group, however using 
azoreductases from Pseudomonas aeruginosa we show that they can rapidly reduce quinones. Azoreductases from the same 
organism are also shown to have distinct substrate specificity profiles allowing them to reduce a wide range of quinones. 
The azoreductase family is also shown to be more extensive than originally thought, due to the large sequence divergence 
amongst its members. As both NAD(P)H quinone oxidoreductases and azoreductases have related reaction mechanisms it is 
proposed that they form an enzyme superfamily. The ubiquitous and diverse nature of azoreductases alongside their broad 
substrate specificity, indicates they play a wide role in cellular survival under adverse conditions. 
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Introduction 

Quinones are best Imown for playing an important role in the 
electron transport chain of most organisms either in the form of 
ubiquinones or menaquinones. The quinones involved in the 
electron transport chain are hydrophobic and lipid soluble. In 
contrast, water soluble quinones are cytotoxic. 

The cytotoxicity of water soluble quinones stems from two 
primary mechanisms [1], firstly their ability to covalently modify 
DNA and proteins, secondly their ability to undergo redox cycling. 

Water soluble quinones are used as defence mechanisms by a 
range of organisms. Water soluble quinones are produced as part 
of the innate immune response of invertebrates via the action of 
pro-phenoloxidase [2] . Plants and fungi use a functionally related 
family of enzymes called polyphenol oxidases (PPOs) to generate 
quinones [3] . The generation of quinones by PPO in plants and 
invertebrates leads to melanisation at the point of infection [4] . 
Some plants also secrete water soluble quinones into the soil to 
inhibit the growth of competing plants e.g. juglone from black 
walnut trees [5]. Bacteria themselves produce quinones as 



antibiotics e.g. doxorubicin from Streptomyces peucetius [6]. Thus it 
is important for bacteria to be able to detoxify these cytotoxic 
compounds. 

One of the primary mechanisms of quinone detoxification is 
their two electron reduction to the more stable and less mutagenic 
quinol [7]. Two electron reduction reactions of quinones are 
carried out by a group of enzymes known as NAD(P)H quinone 
oxidoreductases. Azoreductases are a related but previously 
thought to be distinct group of NAD(P)H dependent flavoenzymes 
that have been identified in a number of species from bacteria [8— 
10], to mammals [11] (Fig. 1). Although azoreductases are 
constitutively expressed in a range of bacteria [10,12], their 
physiological role remains unclear. Azoreductases are able to 
reduce a variety of substrate classes and some have been shown to 
reduce quinones however they have been tested against only a 
limited range of these substrates [8,13]. 

Azoreduction has been shown to proceed via a bi-bi ping pong 
mechanism [14]. An X-ray crystallographic study into the 
mechanism of azoreduction by these enzymes suggested that the 
azo substrate must tautomerise to a hydrazone conformation prior 



PLOS ONE I www.plosone.org 



1 



June 2014 | Volume 9 | Issue 6 | e98551 



Azoreductases and NQOs Belong to the Same Superfamily of Enzymes 



hNQ01 

rNQ01 

hNQ02 

rNQ02 

PA4975 

PA1224 

ecMdaB 

pa2580 

PA1225 

paAzoRI 

efAzoR 

paAzoR2 

paAzoR3 

ecAzoR 

PA 1204 

atNQO 

ecYieF 

pnNQO 

ArsH 

PA2280 

shNQO 

bsAzoR 

rsAzoR 

erNQO 

reNQO 

xaAzoR 

dgFlav 

ecFlav 



I PA0949 

I pcNQO 

I ecWrbA 

I trtiNQO 

Figure 1. Phylogenetic tree illustrating relationships between known and suspected azoreductases from P. aeruginosa PA01. An 

unrooted bootstrap consensus maximum parsimony tree was generated using 500 replicates in Mega 6 [78]. This tree was based upon a sequence 
alignment of 33 proteins using Muscle [79]. The tree was then rooted so that it is consistent with three previously published trees [74,80,81 ]. Enzymes 
in red are from P. aeruginosa, those in green are other bacterial enzymes, those in blue are mammalian enzymes, those in purple are plant enzymes, 
those in pink are fungal enzymes and those in yellow are archeal. PA0949, PA1204, PA1224, PA1225, PA2280, PA2580, and PA4975 are proteins from 
P. aeruginosa. ecAzoR, bsAzoR, efAzoR and rsAzoR are azoreductases from £ coli, Bacillus subtilis, Enterococcus faecalus and Rhodobacter sphaeroides. 
hNQOI hNQ02, rNQOI and rNQ02 are human and rat azoreductases. xaAzoR is a flavin independent azoreductase from X. azovorans. ecMdaB, ecYieF 
and ecWrbA are NAD(P)H quinone oxidoreductases from E. coli. afNQO, pnNQO, tmNQO, pcNQO and atNQO are NAD(P)H quinone oxidoreductases 
from Archaeoglobus fulgidus, Paracoccus denitrificans, Triticum monococcum, Phanerochaete chrysosporium and Arabidopsis thaliana respectively. ArsH 
is an azoreductase from Sinorhizobium meliloti. dgFlav and ecFlav are flavodoxins from Desulfovibrio gigas and E. coli respectively. shNQO, reNQO and 
erNQO are uncharacterised proteins from Staphylococcus haemolyticus, Ralstonia eutropha and Erwinia chrysanthemi. 
doi:1 0.1 371 /journal.pone.0098551 .gOOl 



to reduction [15]. In this meclianism reduction occurs via liydride 
transfer from the flavin group of the enzyme to a quinoneimine 
formed via tautomerisation of the azo bond. This mechanism 
therefore suggests that hydride transfer would also be possible to a 
quinone substrate [16]. Results from studies into reduction of 
quinones by a human azoreductase support this hypothesis by 
showing that that reduction occurs via an obligate two-electron 
reduction [17,18] rather than via one-electron transfer as occurs in 
NADH dehydrogenases". 

P. aeruginosa is known to be a highly versatile pathogen, able to 
infect a range of hosts including invertebrates [19], plants and 
mammals [20]. Three azoreductases, called paAzoRl-3, have 
been characterised in P. aeruginosa [15,21]. Prehminary data 
suggest that paAzoR 1 from P. aeruginosa plays an important role in 
mammalian infection by the bacterium [22]. In light of the 
importance of quinones in the response of plants [23] and 
invertebrates [24] to infection by pathogenic Pseudomonads, 
azoreductases may also play a vital role in infection of other 
organisms via detoxification of their toxic quinones. We have 
therefore carried out a systematic study of the substrate specificity 
profiles of three recombinant azoreductases from P. aeruginosa. W e 
also predict that the azoreductase family is much more extensive 
than originally thought. This data not only improves our 
understanding of the possible role of azoreductases in host/ 
pathogen interactions but, also helps shed light on genomic data. 



Experimental Procedures 

Bioinformatics 

The flow-diagram for the identification of azoreductase 
homologues in P. aeruginosa is provided in Fig. SI. An extensive 
literature review was performed via Google Scholar (search term: 
azoreductase) to identify proteins that have been characterised as 
able to reduce azo compounds. Non-redundant protein sequence 
database search was performed via Blastp [25] to identify enzymes 
related to these azo reducing proteins in the P. aeruginosa PAOl 
genome [26] . Sequences of protein homologues should align well 
to at least 80% of the search sequence. In parallel proteins with a 
similar overall fold to paAzoRl (Pfam families: flavodoxin_2 and 
FMN_red) were identified via the "3D similarity" tool within the 
PDB [27]. Literature searches were then performed on the 
proteins identified to have the same overall fold as paAzoRl, to 
determine which have a characterised quinone reductase activity. 
Blastp searches were then performed using the NAD(P)H quinone 
oxidoreductases identified to search for homologues in the P. 
aeruginosa PAOl genome. Three previously distinct protein families 
(ArsH (named for arsenic resistance), tryptophan repressor binding 
protein A (WrbA) and modulator of drug activity B (MdaB)) were 
identified as putative azoreductases via these searches. The 
original search terms were thus expanded to include these families 
in order to help identify characterised proteins more closely related 
to the putative azoreductases in P. aeruginosa. All sequence 
similarities quoted in the text are based upon pairwise Blastp 
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alignments. Once the azoreductase homologues had been 
identified in P. aeruginosa PAOl, their sequences were used to 
perform searches in other Pseudomonas genomes using Blastp. 

Flavin-independent azoreductases such as AzoB from Xenophilus 
azovorans [28], were excluded from this analysis. This decision was 
made because they are a structurally unrelated (Pfam family: 
NAD_binding_10) family of monomeric enzymes. FMN indepen- 
dent azoreductases use an unrelated mechanism to reduce azo 
substrates that requires them to bind NAD(P)H which directly 
reduces the azo dye bypassing the intermediate steps of FMN 
reduction and oxidation. 

Enzymatic Assays 

All reagents were obtained from Sigma-Aldrich. The pure 
recombinant proteins paAzoRl (PA0785), paAzoR2 (PA1962) and 
paAzoRS (PA3223) were expressed and purified as described 
previously [15,21]. Rates of quinone reduction were obtained by 
monitoring the absorbance at 340 nm for oxidation of either 
NADPH (paAzoRl) or NADH (paAz()R2 and paAzoRS). 
Reactions were carried out in a 100 |J.L reaction volume in UV 
transparent 96-well plates (Greiner) and monitored using a 
Flurostar Omega plate reader (BMG labtech). Reaction mixtures 
were setup as follows: 50 |iM quinone, 500 |lM NAD(P)H and 
between 10 |a,g and 0.1 Hg enzyme in 20 mM Tris-HCl pH 8, 
100 mM NaCl, 5% (v/v) DMSO. Controls were performed in the 
absence of enzyme and all reactions were initiated by addition of 
95 |i,L of enzyme and NAD(P)H to 5 |a,L quinone. Rates were 
determined by fitting the change in OD34() over the first five 
minutes. The quinones tested were as follows (Fig. 2): Coenzyme 
Qj (UQ,1), Phenol blue (Phb), Benzoquinone (Bzq), Juglone (Jug), 
Lawsone (Law), Menadione (Men), Plumbagin (Plu), Adreno- 
chrome (Adr), 1,2-napthoquinone (Onq), 2,3-dichloro-l,4-naptho- 
quinone (Den), 2,5-dichlorobenzoquinone (Deb), 2,3,5,6-tetra- 
chloro-l,4-benzoquinone (Tcq), 2,6-dichloroquinone-4- 
chloromide (Ibc) and anthraquinone-2-sulphonate (AQN). AU 
quinones were prepared as 20 mM stocks dissolved in DMSO. 

Structure Determination 

Solutions of paAzoRl with either UQ_1 or AQ_N were made up 
as foUows: 23 mg.mL"' paAzoRl in 20 mM Tris-HCl pH 8, 

2 mM AQ!SJ/UQ,1 and incubated for one hour at room 
temperature before being washed with 20 mM Tris-HCl pH 8 



to remove DMSO. Crystals grew isomorphously under the same 
conditions described previously [21]. Data for paAzoRl-AQlV 
were collected at the Diamond Light Source (Oxon, UK) on 
beamline 103 using a Quantum ADSC detector. Data for 
paAzoRl-UQ,l were collected at the ESRF (Grenoble, France) 
on beamlin(; ID23.2 with a MAR225 detector. paAzoRl-AQN 
and paAzoRl-UQ,l data processing were carried out in xia2 [29] 
using XDS [30], and SCALA [31]. The structures were solved via 
molecular replacement using Phaser [32], and the structure of 
paAzoRl stripped of its ligands (PDB: 2V9C [21]) as the search 
model. The atomic models were rebuilt in Coot [33] and refined 
using Refmac [34] and PHENIX [35]. Restraints for AQN and 
UQ^l were generated in eLBOW [36]. Model validation was 
performed with MolProbity [3 7] . Data collection and refinement 
statistics are shown in Table 1 and the structures were deposited in 
the protein databank (PDB ID: 4N65 and 4N9Q). 

Results and Discussion 

Quinone Substrate Specificity Profiles 

The substrate specificity profiles of paAzoRl, paAzoR2 and 
paAzoR3 are shown in Figure 3. Each azoreductase has a distinct 
substrate specificity profile. The quinones/quinoneimines tested 
were either benzoquinones or napthoquinones with both ortho 
and para cjuinones (Fig. 2). Interestingly paAzoR3 has a quinone 
reductase activity (Plu -16.8 mM.s '.mg protein ') ~30 fold 
higher than that of paAzoRl (Bzq -0.57 mM.s '.mg protein '). 
Similarly, the highest activity for paAzoR2 is ~ 1 3 fold higher (Bzq 
-7.78 mM.s~\mg protein"') than that of paAzoRl. Rates of 
quinone reduction in all cases were up to two orders of magnitude 
higher than were obser\'ed for reduction of azo substrates as the 
maximum rate observed was paAzoR3 reduction of methyl red 
91.4 |lM.s~'.mg of protein"' [15], suggesting the enzymes are 
better suited for carrying out quinone rather than azo reductions. 
Determination of the Km and Vmax values for quinone reduction 
were hampered by the poor aqueous solubility of the quinones and 
the twin limitations on NAD(P)H concentration to remain within 
the linear range of the plate reader and maintain a greater than 
5:1 molar ratio of NAD(P)H to quinone. We were able to confirm 
that the specific activities measured for both Bzq (paAzoRl and 
paAzoR2) and Plu (paAzoRS) were within the initial linear part of 
the rate curve These data confirm that all three azoreductases are 




Figure 2. Structures of the quinone substrates tested in this manuscript. (A) Bzq, (B) Deb, (C) UQ1, (D) Tco, (E) Phb, (F) Ibc, (G) Horn, (H) Den, 
(I) Men, (J) Law, (K) Plu, (L) Jug, (M) Onq, (N) Ade, (0) AQN, (P) Co7. 
dol:l 0.1 371/journal.pone.0098551 .g002 
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Table 1. Processing and refinement statistics for paAzoRl complexed witli AQN and UQl. 





Structure 


paAzoR1-UQ1 


paAzoRI-AQN 


Space Group 


P3, 2 1 


P3, 2 1 




90, 90, 1 20 


90, 90, 1 20 


0. b, c (A) 


82.40, 82.40, 108.47 


81.89, 81.89, 109.32 


Processing Statistics 


Resolution Range (A) 


43.18-2.0 


43.3-1.8 


Unique reflections^ 


26,575 (2,687) 


40,071 (2,713) 


^merge 


0.053 (0.462) 


0.061 (0.351) 


<I/CT{I)>= 


17.1 (2.1) 


15.5(3.2) 


Completeness %^ 


90.7 (64.0) 


94.9 (66.0) 


Multiplicity'' 


4.7 (2.5) 


5.1 (3.8) 


Refinement Statistics 


Rwork % 


18.51 


16.24 


Rfree % 


23.57 


18.51 


RIVIS bond angle (°) 


0.714 


1.3 


RIVIS bond length (A) 


0.003 


0.011 


Ramachandran Statistics'' 


Preferred region % 


98.5 


98.5 


Allowed region % 


1.5 


1.5 


Outliers % 


0.0 


0.0 



^Numbers in parentheses are for the highest resolution shell. 
"^Ramachandran statistics were calculated using MolProbity [37]. 
doi:1 0.1 371 /journal.pone.0098551 .tool 



NAD(P)H quinone oxidoreductases, however they also highlight 
that the azoreductases are not all equally well adapted for the role. 
The large difiFerences in quinone reduction rates may be partly due 
to the redox potential of the FMN group of the three enzymes. In 
paAzoR2 and paAzoRS the negative charge imparted to the FMN 
after reduction by the NAD(P)H is predicted, based upon 
homology models [15], to be stabilised by interactions with 
His 144 in the case of paAzoR2 and Tyrl45 in the case of 
paAzoR3. In the case of paAzoR2 the interaction with His 144 is 
supported by the recently published native structure of its 
homologue (ppAzoR) from P. putida MET94 (7 1 % identical 
PDB: 4C0W [38]). In contrast in paAzoRl there is no such 
stabilisation as the equivalent residue is Phel50 which has no polar 
group to interact with FMN. 

paAzoRl has a preference for benzoquinones e.g. Bzq (Fig. 3a). 
In contrast paAzoR3 prefers naphthoquinones e.g. Plu (Fig. 3c). 
paAzoR2 reduces both benzoquinones and naphthoquinones 
efficiently (Fig. 3b). One contributing factor to the substrate 
specificity may be structural, paAzoR2 and paAzoR3 are thought 
to have significantly larger active sites than paAzoRl. The 
increased size of the active site of paAzoR2 is supported by the 
native crystal structure of its homologue ppAzoR [38] where the 
heUx bearing Phe60 shifts significantly outwards, away from the 
FMN, thus increasing the volume of the active site. In the case of 
paAzoR3, there are no close homologous structures, however a 
homology model [15] predicts a F60A substitution that would 
significantly increase the size of the active site. A substitution at 
this position was also shown to affect the substrate specificity of 
AzoR from E. coli [39]. The larger size of the active site results in 
paAzoR3 being able to accommodate significantly larger quinone 
groups in an orientation similar to that observed for AQN binding 
to paAzoRl (Fig. 4a). This data shows that NAD(P)H quinone 



oxidoreductases from the same organism have complimentary 
substrate specificity profiles. This had been hinted at by previous 
ad hoc studies of individual enzymes from E. coli (Table 2). This 
substrate specificity data is also consistent with the results from an 
earher study on azoreduction (Fig. S2 [15]). 

Another factor to consider is the redox potential of the 
quinones, as a previous study has shown that there is a linear 
relationship between the redox potential of nitroaromatic com- 
pounds and their rate of reduction by rat NQO 1 [40] . Qiiinones 
have a range of redox potentials, Bzq has a redox potential of + 
0.696 V, while AQN has a redox potential of +0.187 V [41]. The 
relatively high redox potential of Bzq may explain why all three 
enzymes reduce Bzq rapidly while the low redox potential of AQN 
may explain why it is a relatively poor substrate for all three. The 
position as well as presence or absence of electron withdrawing 
and donating groups also changes the redox potential of quinones 
(e.g. Jug has a redox potential of +0.452 V while Law has a redox 
potential of +0.362 V [41]) and may explain some of the large 
differences observed between structurally related quinones e.g. Jug 
and Law in the case of both paAzoR2 and paAzoRS (Fig. 3b and 
3c). 

Co7, was recentiy identified as an inhibitor of human arylamine 
N-acetyltransferase 1 [42]. This data would be of potential 
importance were Co7 to be developed as a therapeutic for breast 
cancer [43]. 

Implications of Quinone Reduction 

Pseudomonads are unusual in being able to colonise a wide 
range of hosts from amoeba to humans. Toxic quinones form an 
integral part of the innate defence of a range of organisms 
including plants, funghi and invertebrates. Several organisms 
produce quinones in response to P.seudomonas infection via the 
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Humans produce melanin primarily via polymerisation of 
DOPA-quinone, in contrast plant PPOs can oxidise a range of 
polyphenols to quinones [47]. We propose the diversity of 
quinones that are produced in turn drives evolutionary pressure 
for bacteria to have a range of NAD(P)H quinone oxidoreductases. 
Plu for example is a potent antibacterial quinone secreted by black 
walnut trees, that is known to have poor potency against P. 
amigmosa [48]. Plu is the preferred substrate of paAzoR2 and 
paAzoRS (Fig. 3). P. fluorescens a close relative to P. aeruginosa that 
has homologues of paAzoRl-3 (Table 3) is insensitive to Jug [49] 
which is rapidly reduced by paAzoR2 (Fig. 3b). 

The hypothesis that the number of azoreductases is related to 
the host range is supported by data showing Pseudomonads that 
infect a small range of hosts have fewer azoreductases than P. 
aemginosa for example P. syringae pv. tomato, a pathogen infecting 
tomato and A.rabidop.sis thaliana, has two azoreductase4ike genes 
(Table 3) in its genome [.50]. A more comprehensive bioinfor- 
matics study of azoreductasc distribution in pathogenic and non- 
pathogenic species, alongside further work with gene deletion 
mutants, will need to be performed in order to test this hypothesis. 

P. aemginosa also produces water soluble melanins with 
benzoquinone functional groups [,^1] which improve the infectiv- 
ity of the bacterium [,52] . During d(;gTadation of phenylalanine 
and tyrosine, Pseudomonads generate the quinol homogentisate 
(Horn [53]). Hom can be readily oxidised to a reactive quinone 
that polymerises to generate melanin [52]. The bacteria needs to 
maintain Hom in its quinol form prior to secretion to prevent 
intracellular polymerisation, hence a requirement for NAD(P)H 
quinone oxidoreductases. All three azoreductases are able to 
reduce a range of l)cnzoquinones (Fig. 3) hence they would also 
reduce the oxidised form of Hom. 




Quinone substrates 



Figure 3. Quinone substrate specificity profiles of paAzoRI (A), 
paAzoR2 (B) and paAzoR3 (C). All rates are normalised relative to 
the maximum rate of NADPH (paAzoRI) or NADH (paAzoR2/paAzoR3) 
oxidation observed for that enzyme. These maximum rates were as 
follows: paAzoRI and paAzoR2 reducing Bzq 0.57 miVl.s~\mg of 
enzyme"' and 7.78 mM.s~\mg of enzyme"' respectively and paAzoR3 
reducing Plu 16.8 mM.s~\mg of enzyme V All rates represent the 
average of three measurements ± standard deviation. 
doi:1 0.1 371/journal.pone.0098551 .gOOB 

oxidation of catechols by PPO. Studies have shown the 
importance of PPO in plant defence against infection by 
pathogenic Pseudomonas syringae [23,44]. Non-pathogenic Pseudomo- 
nas fluorescens has also been shown to induce PPO expression in rice 
plants [45]. The fungus Agaricus hisponis up-regulates PPO 
expression in response to Pseudomonas tolaasii infection [46] and 
in arthropods there is also a significant up-regulation of pro- 
phenoloxidase in response to infection by P. aeruginosa [24] . These 
results indicate quinone production is an important response to 
infection by pseudomonads by many hosts. 



Structures of paAzoRI Binding UQ1 and AQN 

In order to help understand the mechanism of quinone 
reduction by azoreductases the structures of paAzoRl bound to 
two quinone substrates (UQ_1 and AQjNf) have been solved via X- 
ray crystallography. The structure of paAzoRl bound to AQjNF 
was solved at a resolution of 1.8 A with final Rj„,, and Ru,„rt values 
of 18.5'% and 16.2% respectively, while UQ_1 was solved at 2.0 A 
with final i?/^^, and Ri„,„k values of 23.6% and 18.5% respectively 
(Table 1). No significant conformational changes were observed 
upon binding of either UQ.1 or AQN to paAzoRl (RMSD 0.22 A 
compared to PDB: 2V9C). In both structures, the quinone is only 
observed binding to one of the two active sites, although positive 
difference densit\' is observed in the alternate active sites there was 
insufficient to unambiguously buUd in the model. During the 
refinement of the structure of paAzoRl hound to UQ_1, the 
occupancy of the UQ_1 molecule was reduced to 0.5 suggesting 
that VQl was not present in £iU paAzoRl molecules. 

In the structure of paAzoRl binding to AQN, AQN makes 
hydrophobic contacts with Phe60, PhelOO, Phel20, Tyrl31 and 
Phel73 (Fig. 4a). AQjN also forms %-% stacking interactions with 
the central ring of FMN while the sulphonate group points into a 
water filled channel. Neither quinone oxygen is in a position to 
form a hydrogen bond with Asn99 upon reduction as has been 
previously observed [54] most probably due to the bulky nature of 
AQN. Like AQN, UQl forms n-n stacking interactions with FMN 
and makes hydrophobic contacts with PhelOO, Phel20, Tyrl31 
and Phel73 (Fig. 4b) but it also hydrogen bonds with Tyrl31 via 
the quinone oxygen and forms extra hydrophobic contacts with 
PhelSl. 

The binding orientations of UQl and AQN differ significantly, 
with UQ^l binding parallel to the ring system of FMN while AQ]S[ 
binds at almost a 90° angle to FMN (Fig. 4). These differing 
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Table 2. Structural and functional information on all azoreductase homologues. 



Protein 


Structure* 


Enzymatic activities 










Azoreductase 


NAD(P)H quinone oxidoreductase 


Nitroreductase 


paAzoRl 


2V9C 


[21] 


This paper 


[54] 


paAzoR2 


ND 


[15] 


This paper 


ND 


paAzoRB 


ND 


[15] 


This paper 


ND 


PA0949 


IZWL 


ND 


[75] 


ND 


PA1204 


1X77 


ND 


ND 


ND 


ecAzoR 


2Z98 


[8] 


[8] 


[83] 


ecMdaB 


2B3D 


ND 


[66] 


ND 


ecWrbA 


2R96 


ND 


[74] 


ND 


ecYieF 


3SVL 


ND 


[84] 


ND 


ArsH 


2Q62 


[58] 


ND 


ND 


rsAzoR 


ND 


[85] 


[86] 


ND 


bsAzoR 


INNI 


[13] 


[13] 


ND 


hNQOl 


IQBG 


[11] 


[11] 


[11] 


hNQ02 


2QR2 


[11] 


[11] 


[11] 


pnNQO 


3U7R 


ND 


[87] 


ND 


atNQO 


ND 


ND 


[70] 


ND 


afNQO 


ND 


ND 


[74] 


ND 


^PDB codes for structures, ND not determined. Where particular 
has been reported. 


enzymatic activities have been observed a 


reference is provided, otherwise no reference to the activity 
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binding orientations correspond to those observed when methyl 
red binds to either the wild type enzyme [21] in the case of AQN 
and the Y131F mutant [14] in the case of UQl. The binding 
orientation of UQ^l is also consistent with the binding of both 
balsalazide and nitrofurazone to the wild type enzyme [15,54]. 
The binding orientation of AQN to paAzoR 1 is consistent with 
bindmg of menadione to hNQP2 (PDB 2QR2). A recent X-ray 
crystal structure shows the binding of AQN to the ppAzoR (PDB: 
4C0X [38]). AQN is shown to bind in an orientation similar to 
that observed for UQ_1 bound to paAzoRl, with the sulphonate 
positioned in solvent outside the pocket. In contrast to the 



structure with paAzoRl, AQN is accommodated inside the active 
site of ppAzoR via a number of significant side chain conforma- 
tional shifts. Due to the altered binding conformation of AQN, 
compared to the paAzoRl structure, a water bridged hydrogen 
bond can form between AQN and Asn97 (Asn99 in paAzoRl). 
This water is similar to that seen in the structure of paAzoRl 
bound to nitrofurazone [54] and is likely to be involved in the 
mechanism of reduction. 

In the case of UQ_1 a 8""" methoxy substituted carbon of the 
quinone ring is 3.5 A from the N5 of FMN making it an ideal 
candidate for hydride transfer which corroborates the previously 




Figure 4. Binding of paAzoRI to AQN and UQ1. Binding of paAzoRI to AQN (A) and UQ1 (B). AQN is shown with pink carbon atoms, UQ1 with 
grey carbon atoms and FMN with yellow carbon atoms. The green mesh is unbiased positive difference density (F^-F^ map) contoured at 3a while the 
blue mesh is the refined 2Fo-Fc map at 1 a. A red dashed line signifies a hydrogen bond. Residues interacting with the ligands are labelled. (A) is based 
upon PDB 4N65, while (B) is based upon PDB 4N9Q. 
doi:10.1371/journal.pone.0098551.g004 
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Table 3. Distribution of azoreductase-lil<e enzymes across Pseudomonas species. 



Gene 


P. fluorescent 


P. putid^ 


P. stutzerf 


P. denitrifican^ 


P. aureofacient 


P. syringae pv. tomato' 


pa0785i' 


78 






68 


77 




pa0949 


84 


83 


86 


92 


81 


81 


pa1204 - - . - 


pal 224 




45 


76 




75 




pal 225 


68 








72 




pal 962^ 


67 


71 




76 


67 


63 


pa2280 


83 


87 & 80 


80 




82 




pa2580 


79 


86 


76 


81 


75 




pa3223' 


69 


67 




74 


66 




pa4975 




78 




82 







All values are given as percentage identity to the protein from P. aeruginosa PA01. Homologues are defined as having >45% sequence identity based upon Blastp 
alignments. The strains whose genomes searched are as follows: 

"pfo-i, 

"^2440, 
'ATCC 17588, 
"ATCC 13867, 
'30-84 and 'DC3000. 

The genes encoding paAzoRI, ^paAzoR2'^ and paAzoR3' are marked. 
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proposed mechanism of quinone reduction based upon compar- 
ison witli azo reduction [16]. In tlie case of AQN one would 
expect the 5"^ carbon, which accepts the hydride, to be the carbon 
attached to the negatively charged sulphate. However the sulphate 
of AQN is not easily accommodated within the primarily 
hydrophobic active site of paAzoR 1 . As a result an unmodified 
carbon atom (Fig. S3) accepts the proton in a less favourable 
reaction possibly contributing to the lower rate of reduction 
observed of AQN by all azoreductases compared to other 
quinones (Fig. 3). 

Identification of the Extended Azoreductase Family of 
Enzymes 

The availability of the sequence of the P. aeruginosa PAOl 
genome [26] allowed initial identification of paAzoRl, paAzoR2 
and paAzoR3 via sequence homology to E. coli AzoR [2 1] using a 
sequence identity cut-off of 30%. Literature reports show that 
proteins with azoreductase and NAD(P)H quinone oxidoreductase 
activities have highly diverse amino acid sequences (Fig. 1). Due to 
their conserved structure (Fig. 5) and mechanism of action [16] we 
propose that the azoreductases and NAD(P)H quinone oxidore- 
ductases form an enzyme superfamily. This sequence diversity 
stems from the fact that all reactions occur on the flavin cofactor 
(either FMN or FAD) which is anchored to the protein via a large 
number of sequence independent interactions with the protein 
backbone [15]. Due to the sequence diversity several putative 
azoreductases were missed during the original bioinformatics scan. 
The following genes are proposed to encode azoreductases: 
pa0949, pal204, pal224, pal225, pa2280, pa2580 and pa4975, 
based upon sequence homology to proteins of known structure and 
function. The phylogenetic tree that illustrates their evolutionary 
relationship to each other and to characterised enzymes from 
other species is shown in Figure 1, the available structural and 
functional data for each characterised enzyme in the tree is 
displayed in Table 2. 

Homologues of known azoreductases. The homologues of 
four putative azoreductase genes from P. aeruginosa have been 
characterised as being able to reduce azo dyes. The first of these is 



the protein encoded by pa22S0 which is a member of the ArsH 
enzyme family. ArsH was initially identified as part of the arsenic 
resistance operon of Tersiniae enterocolitica [55] although gene 
deletion does not affect resistance to arsenic. A functional arsenic 
resistance operon was identified in P. aeruginosa PAOl [56] 
consisting of Ai-sR, ArsB and Ai-sC (PA2277-PA2279) and 
pa2280 is predicted by DOOR [57] to be part of the same operon 
as these genes. ArsH from both Sinorhizobium meliloti (82% similarity 
[58]) and Synechocystis sp PCC 6803 (78% similarity [59]) have been 
shown to reduce azo dyes. Synechocystis sp PCC 6803 ArsH also 
reduces quinones [59]. The structure of ArsH from S. meliloti [58] 
has the same overall core fold as paAzoRl (RMSD = 2.2 A — 
Fig. 5a). The major difference in structure between paAzoRl and 
the ArsH family is a C-terminal extension of around 30 amino 
acids that is absent from most other members of the azoreductase 
family. The function of this C-terminal extension is unclear. 

The second putative azoreductase in this class is encoded by 
pa4975. The gene product of pa4975 is part of the arginine 
transaminase operon of P. aeruginosa [60] although its function 
within this operon remains unclear. PA4975 is unusual in bacteria 
in that it shows significant sequence similarity to the human 
enzyme NQ02 (61% similarity). hNQP2 regulates proteasomal 
degradation of a range of proteins [61], and also regulates TNFa 
induced apoptosis [62]. hNQP2 can reduce both azo and quinone 
substrates [1 1]. The structure of hNQP2 has been solved and has 
the same overall fold as paAzoRl (RMSD = 1.5 A - Fig. 5a). 

The third and fourth azoreductase-like enzymes are encoded by 
pal 224 and pal 225 which show sequence homology to the human 
azoreductases hNQOl (54% similarity) and hNQP2 (48% 
similarity) respectively. hNQOl also known as DT-diaphorase, 
regulates proteasomal degradation of various proteins including 
transcription factors such as p53 and C/EBPot [63]. Although 
adjacent on the chromosome the proteins are encoded on opposite 
strands and hence are unlikely to be co-regulated and little is 
known of the function of either gene. hNQOl, like hNQP2, is 
both an azo and quinone reductase [11]. hNQOl shares the same 
overall fold as paAzoRl (RMSD 2.1 A [64]) but like ArsH has a 
long C-terminal extension. 
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Figure 5. Structures of azoreductase homologues overlaid onto paAzoRL (A) Structures of ArsH (PDB: 2Q62 - lilac), ecMdaB (PDB: 2B3D - 
coral) and hNQ02 (PDB: 2QR2 - brown) overlaid onto paAzoRl. (B) Structures of PA0949 (PDB: IZWL - red) and PA1204 (PDB: 1 x77 - green) overlaid 
onto paAzoRI. In all cases, one monomer of the dimer is shown. In both cases, paAzoRI is in blue (PDB: 2\/9C). All alignments were carried out via 
secondary structure matching and images were generated in CCP4MG [82]. 
doi:10.1371/journal.pone.0098551.g005 



Homologues of known NAD(P)H quinone 
oxidoreductases. Close homologues of several other P. aerugi- 
nosa genes have an NAD(P)H quinone oxidoreductase activity and 
the same overall fold as paAzoRl but have yet to be tested for 
azoreductase activity. The protein encoded by pa2580 shares 
significant homology with a group of enzymes known as MdaB. 
The MdaB family was initially identified in E. coli due to its 
overexpression causing increased doxorubicin resistance [65]. The 
structure of i?. coli MdaB (75% similarity) has been solved and has 
the same overaU fold as paAzoRl (RMSD = 2.2 A - Fig. 5a). 
MdaB from both E. coli [66] and Helicobacter pylori (7 1 % similarity 
[67]) have NAD(P)H quinone oxidoreductase activities. 

The protein encoded by pal 204 is part of the BexR regulon in 
P. a£ruginosa [68] and is part of an operon consisting ot pal 202- 
pal205 that is strongly up-regulated in BexR gene {pa2432) 
deletion mutants. The structure of PA1204 has been solved and 
was shown to have the same overall fold as paAzoRl 
(RMSD = 2.3 A -Fig. 5b [69]). The closest characterised homo- 
logue of PA1204 is an enzyme from A. thaliana that shows 
NAD(P)H quinone oxidoreductase activity (69% similarity [70]). 

The protein encoded by pa0949 is part of the WrbA family, 
named for their strong interaction with tryptophan repressor [71]. 
In P. aeruginosa as in many bacteria, expression of PA0949 is 
controlled by the PhoP/PhoQ_ two component regulatory system 
[72]. The structure of paWrbA has been solved and shows the 
same overall fold as paAzoRl (RMSD = 2.3 A - Fig. 5b [73]). The 
E. co/j homologue of WrbA (56% similarity) has NAD(P)H quinone 
oxidoreductase activity [74]. Although phylogenetic analysis does 
not support a common ancestor of the azoreductase and WrbA 
enzyme families (Fig. 1), we believe the structural and enzymatic 
data is sufficient for it to be included as a likely candidate to be an 
azoreductase. 

Recent data have shown that PA0949, PA1204, PA2580 and 
PA4975, have NAD(P)H quinone oxidoreductase activities [75]. 
This data supports their inclusion in this list of possible 
azoreductases. 

Distribution of azoreductase like genes in Pseudomonas 
species. It is interesting to note that not all Pseudomonads have 
the same complement of azoreductases that are observed in P. 
aeruginosa PAOl (Table 3). A range of genomes from both 
fluorescent and non-fluorescent Pseudomonas species were searched 
for proteins that were close homologues (>45% sequence identity) 
of those found in P. aeruginosa PAO 1 . In the six genomes that were 



selected only one putative azoreductase (paWrbA or pa0949) was 
present in all. WrbA homologues are also found in a range of other 
bacteria including E. coli (Table 2) however WrbA is not essential 
for survival [74]. In contrast, close homologues of PA1204 were 
absent from all six genomes. Interestingly the closest homologue of 
PA1204 in P. Jluorescens pfD-1 (Pfl01_2 192-38% identical) does 
have homologues in P. putida (69% identical), P. denitrificans (74% 
identical), P. aureofaciens (86% identical) and P. syringae pv tomato 
(75% identical). 

Conclusions 

Azoreductases are a varied family of enzymes that have been 
identified in almost all species (with the exception of viruses). Their 
ability to reduce a wide variety of both endogenous and exogenous 
compounds has complicated the identification of their physiolog- 
ical substrate and hence their function. The quinones identified in 
this study appear to be the best substrates for the enzymes yet 
identified, indicating that they may act as the physiological 
substrate for these enzymes. Enzymes with a flavodoxin-like fold 
are frequently observed to have both azoreductase and NAD(P)H 
quinone oxidoreductase activities (Table 2) therefore we suggest 
that these activities are indistinguishable in these enzymes. 

Water soluble quinones are an important part of the defence 
mechanism of many organisms against bacteria including P. 
aeruginosa. In order to detoxify these quinones P. aeruginosa must 
reduce the quinone to the more stable quinol. Plants in particular 
are able to oxidise a range of catechols to quinones via their PPO 
enzymes, and as a result Pseudomonads need to be able to reduce 
a range of different quinones in order to survive in plant systems. 
We have shown the azoreductase family are very active NAD(P)H 
quinone oxidoreductases with broad specificity profiles and as a 
result they may play an important role in detoxilying the quinones 
secreted by many of the host organisms of P. aeruginosa. Infection of 
model plants with azoreductase gene deletion mutants in P. 
aeruginosa could be studied to improve our understanding of the 
role of azoreductases in plant infection. 

In addition to their enzymatic activities a number of 
azoreductase homologues have been shown to alter the action of 
proteins via mechanisms independent of their catalytic function. 
The human azoreductases liNQOl and hNQP2 regulate 
ubiquitin-independent proteasomal degradation of various pro- 
teins via direct protein-protein interactions [61,63], while in E. coli 
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the NAD(P)H quinone oxidoreductase KefF [76] interacts with the 
potassium channel KefC and regulates its ion flux [77]. The 
importance of these protein-protein interactions in the physiolog- 
ical function of azoreductases requires further detailed study. 

Supporting Information 

Figure SI Flowchart outlining the bioinformatics pro- 
cedure applied to identifying members of the azoreduc- 
tase family. 

(PDF) 

Figure S2 Azo compound substrate specificity profiles 
of paAzoRl (A), paAzoR2 (B) and paAzoR3 (C). Abbrevi- 
ations for all substrates are as follows; balsa - balsalazide, pabsa - 
p-aminoazobenzene-4'-sulfonate, sulfa - sulfasalazine, olsa - 
olsalazine, metred - methyl red, tropa- tropaeolin O, amar - 
amaranth, pone S/BS - ponceau S/BS,. All rate.s are normalised 
relative to the maximum rate of reduction of Balsalazide 
(paAzoRl -29.7 (iM.s"'.mg"'), ponceau BS (paAzoR2 - 
26.3 nM.s~\mg"') or methyl red (paAzoRS 
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